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k> ! Abstract 

j^ ■ The KdV and modified KdV integrable hierarchies are shown to be different 

descriptions of the same 2D gravitational system - open-closed string theory. Non- 
perturbative solutions of the multi-critical unitary matrix models map to non- 
singular solutions of the 'renormalisation group' equation for the string susceptibil- 
ity, [P, Q] = Q . We also demonstrate that the large N solutions of unitary matrix 
integrals in external fields, studied by Gross and Newman, equal the non-singular 
pure closed-string solutions of [P, Q] = Q . 
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Introduction. 

In this letter we sliall stiow tiow double-scaling limits of unitary matrix models 
can be interpreted in terms of 2D quantum gravity and its worldsheet expansions. 
The modified KdV and KdV integrable hierarchies are shown to be different descrip- 
tions of the same string theory. Specifically, using the Miura transformation between 
these hierarchies, and the inverse transformation which we explicitly construct, we 
map the solutions of the double-scaling limit of symmetric unitary matrix models 
with C flavours of massless 'quark' [1] [2] to non-perturbative solutions of open- 
closed string theory in the (2, 2m — 1) minimal model backgrounds. The open string 
coupling r, measured in units of the closed string coupling z/, satisfles C = 1/2 ± F 
and the boundary cosmological constant is zero. The map implies that the two 
physical r-functions of the mKdV/KdV hierarchy [3] equal square roots of partition 
functions of open-closed strings in general. As a simple corollary the continuum 
limit 1/A^-expansion of the original unitary matrix model partition function [1] is 
seen to be nothing other than world-sheets with even numbers of boundaries and 
r^ = 1/4. The general form of the Virasoro constraints on the r-function are 
explained. Correlators of local operators can be evaluated in terms of the flows 
of the hierarchy in the usual way. 

Our analysis also leads to a proof of some transformation equations between 
solutions of the string equations with open string couplings differing by an integer, 
generalising a relation proved by Lukashevich, Fokas and Ablowitz[4] for solutions of 
Painleve II. Although the Miura transformation from the mKdV hierarchy applies 
to very general solutions, of particular interest are the unique real non-singular 
solutions furnished by the unitary matrix models themselves [5] [6]. These map 
to open string generalisations of the non-singular solutions of the KdV hierarchy^ 
found in [8] [9] . Remarkably, the pure closed string partition functions uncovered 
in those papers are almost surely the solutions determined by Gross and Newman 
[10] [11] for the simplest unitary matrix model in an external fleld and its formal 
multicritical generalisations. Damning evidence of this fact will be presented in a 
later section. The main results are expanded upon in the following section, but flrst 



^ For the case of coupling to the general (p, q) minimal models see also ref. [7] 
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we conclude this introduction by recalling some salient points about open-closed 
string theory and matrix models. 

Solutions u{zj tk) of the KdV hierarchy^ dt^u = dzTZk+i [u] must satisfy a couple 
of physical requirements if they are to represent the string susceptibility of a 2D 
gravity theory. Firstly as a function of the dimensionful arguments {tk, z} they must 
satisfy the renormalisation group (RG) equation for invariance under a change of 
scale [8] , given by [P, Q] = Q in terms of pseudo-differential operators [9] . After 
some simple manipulation one can rewrite it [8] : 

^t7^2--7^7^" + ^(7^V = ^2 (i) 

2 4 

where ' = dz, TZ is defined in terms of Gel'fand-Dikii differential polynomials 7^fc[w] 
as 7^ = J2k=i{k + l/2)tk'R.k — z, and F^ arises as an integration constant. Secondly u 
must have an asymptotic expansion in ^ as ^ ^ oo which matches that obtained from 
the perturbative genus expansion of 2D gravity coupled to some matter theory. For 
the (flows between the) (2, Iva — 1) minimal models coupled to gravity this means 
that for closed strings u must satisfy 7^ = to all orders in the asymptotic expansion. 
Eqn.(l) has such solutions provided that F = 0. Indeed, either 71 = exactly [12] 
yielding solutions which are complex [13] [14] or have singularities somewhere on the 
flow [15], or one may choose the solution with unique real non-singular flow studied 
in refs.[8][9][16] (referred to as DJM solution hereafter). F plays the role of an open 
string coupling as follows by comparing eqn.(l) with the string equation found by 
Rostov [17], obtained by performing the double-scaling limit of matrix models with 
a term -^yTr log (1 -|- (j)"^) in the potential to generate random holes [18]; each surface 
is weighted by r#hoies -^^ ^^^^ Rostov's solutions were restricted to those for which 

7^ = -2Fdiag{(-a2 + «)-!} (2) 

these obeying 7^ = in the F — i> limit. Hence eqn.(l) is more general, encompassing 
the open string generalisations of the DJM closed string solutions, the previous 
relation (2) then holding as a matching condition in the sense of an asymptotic 
expansion as z ^ oo, but being violated at the non-perturbative level. Note that 



'^ z is essentially the 'space' variable in the RdV equations, the tk with k > 1 
being generalised time variables. 
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the renormalised boundary cosmological constant p which one may assign to the 
holes, determined by the position of the branch point of the logarithm, has been set 
to zero in (1) since this is the case which will interest us in what follows. We shall 
discuss p 7^ briefly later, in the context of the Virasoro constraints. 

Introducing the closed string coupling z/ (renormalised 1/A^) into the string 
equations (1) and (2) by the rescaling tk -^ tkj^, z -^ zjv, the asymptotic solution 
is a series in v which at the mth critical point takes the form: 



,1/m 



^2g+7ip/i 



g,h=0 

Agh is determined in (2) once Aqq has been flxed, and determined in (1) once Aqq 
and the sign of Aqi have been flxed. For the sake of clarity let us rescale such that 
{m + l/2)tm'T^m = W^ + ■ ■ "j then in the case z — > +C)0 the sphere term is required 
to be Aqq = 1. Using (2) then flxes the disc term to Aqi = — l/rri. We shall also 
need later the real z ^ —oo expansion particular to the DJM solution, obtained by 
taking the new possibility allowed by eqn.(l), Aqq = 0: 

-.2 / -I X oo mr 2mrp2s 

^ \ ^^ r=Os=0 ^ 

Substituting (4) in (1), the torus and cylinder terms are readily seen to be as above 
with Bqq = 1. The general form then follows from uniqueness together with the 
symmetry of eqn.(l) (after introducing v) under z/^ ^ J'^^, u ^ ju where j is any 
mth root of one. 

String Equations and the Miura Map. 

Solutions v{z,tk) of the modifled KdV hierarchy are deflned by the flow 
equations dt,.v = ^dzSk[v] where Sk = ^IZ'fJyv'^ + v'] — vTZk[v'^ + v']. Here v"^ plays 
the role of string susceptibility, and its analogous RG equation is just the differential 
of the unitary matrix model string equation [8] given by the flowing version of the 
Painleve II hierarchy: 

OO -1 

E (^ + o^tkSk[v{z)] + zv{z) = C (5) 

k=i ^ 



(C is again an integration constant). This equation with C = was first found by 
Periwal and Shevitz [1] by taking the double-scahng hmit of a unitary matrix model 
with partition function: 

Z= I r't/exp(A^Tr ^ gp{U + U^Y) (6) 

The eigenvalues of U lie on the unit circle and the critical point occurs when the 
support of their large A^ density function just covers 27r. The two ends of this 
support meet to produce 2m zeros when the potential is tuned to the mth critical 
point. The significance of C 7^ was later understood by Minahan [2] who derived 
eqn.(5) by adding a term ±^Tr log |1 + t/| to the potential, which may be obtained 
by integrating out C fiavours of bosonic or fermionic 'quarks'. 

We will now display a one to one map between solutions of (5) with C = 1/2 ± F 
and solutions of (1) . This map is the Miura transformation u = v"^ + v' familiar 
in the theory of the KdV and mKdV equations^. Furthermore we will define a 
one to one transformation between solutions vc{z) of (5) and solutions of (5) with 
unit change in the boundary coupling: vc±i{z). Equivalently between ur{z) and 
u\Y'±i\{z). To avoid confusion we need to state carefully our conventions. The 
maps apply to general solutions ur of (1) and vc of (5), not only those obeying the 
physical boundary conditions discussed below eqn.(2) and below (13). The value 
of the boundary coupling is used to label any one of the set of solutions with that 
coupling, and different labels may refer to different functions of the two arguments 
e.g. setting v±{z, C ± 1) = vc±iiz) we have v+{z, A) 7^ V-{z, A) in general. Finally 
the label F for solutions ur of (1) will always be taken to be no n- negative; the actual 
sign of the open string coupling in the physical solutions is fixed by reference to (2) 
and will be discussed later. 

Let u{z) be a solution to (1) and define 

X±[u,v] = ^n[u]TT-v{z)n[u] . (7) 

Then X±[u, v] = defines a function v{z) (possibly with singularities): 

TZu] ^ ' 



^ The authors of ref.[3] used the Miura map in the case 7^ = (F = 0) to obtain 
the Painleve II hierarchy, but incorrectly dropped the constant C — 1/2. 
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Note that one can define v for tfie F ^ limit of Rostov's solutions by using (2) to 
cancel F first. Now 

= xl ± 2rx± - n[u]x'^ = {v^ + v')n'^[u]- -n[u]n"[u] + hn'[u])'^ -r^ . (9) 

By comparison with (1) we deduce that the inverse transformation is u = v'^ + v' , 
which on substituting in (8) and rearranging gives eqn.(5) with C = 1/2 ± F. On the 
other hand given a solution v of eqn.(5) with C = 1/2 ± F, and defining a function 
u = v'^ + v' we have X±[u, v] = hj rearrangement, and hence (8) is the inverse 
transformation and u satisfies (9) which is (1). 

Thus we have constructed a one to one correspondence between solutions vc 
of (5) with C = 1/2 ± F and u-p of (1) . Since this is true for both choices of C 
it follows from the Miura map that there is a one to one correspondence between 
solutions Vc and solutions vi-c such that 

Vc + v'c = vl_c + v[_c 
and by using (8) 



2C-1 

nvh+^'c] 



Vi_c = Vc+ ^r. 2 , ../ 1 ' (10) 



It is easy to see that this 'fiip' transformation changes the sign on boundary 
conditions for v given by f ^ -|-^i/2m ^^ ^ ^ ^^ q,^ ^^le other hand (5) is odd 
in V [1] . One can see this for example by induction, using the linear recurrence 
relation in TZrn to prove 

1 .„, 1 



r>' / O Olll I oil i Cl 



'm+l ~ fJ <^m+l 



4 "^ 4 



m 



and noting that since the TZm [u] contains no -u-independent terms, Sm can contain 
no terms linear in v. Thus, given a solution vc we can construct a solution v-c by 

V-C = -Vc ■ (11) 

Combining this with (10) we obtain the following one to one correspondence between 
solutions which in particular preserves the above boundary conditions: 

2C±1 
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By using the Miura map we obtain similarly a one to one correspondence between 
solutions ur{z) and Uy{z) when F ± F is an integer. It is clear that there are 
many properties that can be deduced about the solutions of (1) and (5) from these 
equations. In particular, generalising ref.[4], if one starts with vo{z) = one can 
generate a series of rational solutions to (5) and (1) with F — 1/2 or C an integer, 
and finitely many tk ^0. 

Returning to 'physics', if we now re-introduce u into (1) and (5) again via tk — ^ 
tfc/z/, z -^ z/u, one finds at the mth critical point that, when it has an asymptotic 
expansion, v is of the form 

oo ^2p+q(jq 

V^±Z "^ Y. '^P9 (2+l/m)(p+g/2) " ^^^) 

p,q=0 

With tm set as below (3), |Coo| = 1 or Coo = 0. In fact the asymptotics v -^ j-^i/Sm 
{v —^0) as z ^ +00 {z -^ — oo) are the ones required by unitary matrix models[l] [2]. 
If C > the sign choice corresponds to C bosonic (fermionic) quarks [2]. (C < 
can be fiipped to these by (11) and in the original C = case[l] the sign is clearly 
irrelevant). For either sign, a simple calculation reveals that Cqi = — 1/(2?ti) as 
z -^ +00, while the leading asymptotic is z; = Cu/z as z ^ — oo. Using the Muira 
map {u = v"^ + udzv) on these low orders of perturbation theory one readily finds 
that the asymptotics map onto those of the DJM solutions as in eqns.(3)(4). Indeed 
this also determines the sign of the open string coupling to be F = ±(C — 1/2), 
again the sign choice being that of the z -^ +oo f-asymptotic. This may be used 
to determine the effect on F of the maps in (10)-(12); the effect in all cases is to 
increase or decrease F by 1, thus generating one to one maps between solutions with 
open string couplings differing by integers. 

Of particular interest are the known non-singular solutions. When C = the 
unitary matrix model (6) provides a unique non-singular fiow between the critical 
points [5] [6] having the plus-branch boundary conditions at the mth critical point. 
Trivially by (11) there are such solutions for the minus branch also. Thus these 
map under Muira onto DJM-type solutions with non-singular fiows, for open string 
couplings F = 1/2 and F = —1/2. Incidentally for z -^ — oo (4) implies that u is 
zero to all orders in the asymptotic expansion, consistent with the f-solutions since 
they also have this behaviour. On the other hand the unique non-singular flows 
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of the DJM solutions [16] for F = map onto the fermion branch^ tJ-solutions with 
C = 1/2. If these are non-singular, which we expect, then they are also the unique 
solutions with this property. 

Generally, there is to our knowledge no thorough study of the existence and 
uniqueness of non-singular solutions of the Painleve II hierarchy with C 7^ in the 
literature, nor for eqn.(l) when F 7^ 0. Employing the techniques used for the C = 
case [5] [19] has not yet borne fruit for us. Minahan argued [2] that there should 
be a unique non-singular real solution for C > fermionic quarks. These would 
map smoothly to real non-singular D JM-type solutions of open-closed string theory 
with F < 1/2. Due to (11) there would also be a C < real pole- free solution with 
positive-branch asymptotic, giving F < —1/2. On the other hand similar arguments 
suggest that there may be unique real non-singular -u-solutions for F < 0. 

This two-fold branching when F 7^ is also reflected in the r-functions. As 
noted in ref. [3] , there are for the joint KdV and mKdV system related by the Miura 
map, two tau-functions tq and ri arising from the two basic representations of the 
Kac-Moody algebra Al . They satisfy 

v'^ ± ud^v = -2u'^dl log To,i (14) 

respectively and are square roots of partition functions of open-closed string theory 
in the (2, 2m — 1) minimal model backgrounds. The unitary matrix partition 
function itself (6) is the product tqTi and so using the fact that tq and ri have open 
string worldsheet expansions with F = ±1/2 respectively (C = 0), the continuum 
worldsheet expansion of (6) has only even numbers of holes and F^ = 1/4. Although 
it might seem odd that this picture of the continuum limit is not all that similar 
to that of the 1/A^ expansion before double-scaling, one must remember that some 
other matrix model with 'two-cut' behaviour of the eigenvalue density would have 
sufficed equally [15] . Rather it should caution us that the universal surfaces 
creamed-off by the double-scaling limit need bear little resemblance to the finite ones 
we see in the Feynman diagrams. It has been suggested that mKdV describes dense 
polymer trees [20]. Though such polymers have a 'dual' description as dense self- 
avoiding loops which is reminiscent of the holes created by open strings, especially 



^ The sign follows from (8) and by inspection of the m — 1 numerical solution. 
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in view of the fact that there is no boundary cosmological constant to suppress their 
length, we have been unable to make a precise identification. Moreover the authors 
of ref.[20] were considering the full 2nd mKP hierarchy, including the 'even' time 
variables. We have been working here with the reduced case obtained by setting 
these to zero and using only the 'odd' time variables tk- 

Before closing this section let us mention the Virasoro constraints of open- 
closed string theory on the tau- functions tq, ti. These are the usual ones of closed 
string theory except that L_i is apparently missing and Lq contains an eigenvalue 
fx (as well as the 1/16) i.e. (L^ — Hidn,o)Ti = where i = 0, 1 and n > [3]. Indeed, 
substituting v = dz log(ri/ro) into the mKdV fiows and using eqn.(5), one readily 
deduces^ that /Ui — /Uq = C/2. Note also that the Lq equation is the RG equation 
for the T function and is equivalent to (1) for solutions of the KdV hierarchy. As 
remarked previously, for open strings we are free to introduce a constant p which 
weights holes of length / in the worldsheet by e~^'. For simplicity the discussion 
which follows is understood to be perturbative^ in u. We can generate r(p) from 
r(0) by using the fact that a boundary cosmological constant couples to L_i [23]; 
L_ir(/9) = dT{p)/dp. This equation may be confirmed by a straightforward matrix- 
model calculation which determines the r.h.s. as 

^ = 2rudzdmg{{-u^dl + u + p)-'} . (15) 

op 

It is then equivalent to 7^ = — 2rz/diag{(— z/^9f + u + p)~^}, which is always true to 
all orders in perturbation theory. The general solution is r(p) = e''^~^r(0). This 
filters through to the Virasoro constraints as 

Lnrip) - eP'^-^i^i-p adL_i)^ ■ L, + -^-—{-p adL_i)"+i ■ L,}r(0) 
n! [n + iy. 

where adL_i • L^ = [L-i, Lm], which on rearranging gives 



LnT{p) 



n+1 9 



Lr,-{n + l)pp^-p^^' — 



r(p) = Vn>-1 . (16) 



^ We thank T.HoUowood for bringing this to our attention [21]. 

^ Amongst other things this obviates the need to burden the reader with certain 
non-perturbative effects on the boundary cosmological constant which arise for DJM 
solution (see [7] [22] for details). 
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It is easy to confirm that these L^ : n > — 1 satisfy the centreless Virasoro algebra. 
We do not wish to dweU upon the details here but merely emphasise the obvious 
but important fact that 

dr 
r^O^ — ^0 for any p (17) 

op 

Thus when p = 0, as we have been considering previously, the constraints at n > 
are the previous ones, while n = — 1 remains there but modified so long as we are 
dealing with open strings. Therefore L_i is not missing but rather is neutralised 
as a constraint by the presence of an extra variable other than the tk- 

Unitary Matrix Model in an External Field. 

In this final section we describe the equivalence of the DJM closed string 
solution (r = 0) to eqn.(l) and that of the N x N unitary matrix model in an 
external field and its formal multi-critical generalisations considered by Gross and 
Newman[10][ll]. This involves the following integral 

Z= I VUex.p{NTY[UA^ + U^ A]) (18) 

which can arise in the mean field approach to lattice gauge theory. A is a complex 
'background' matrix and some Dyson- Schwinger equations for Z were used in [10] 
to solve the model in principle to all orders in 1/A^, the model depending only 
on the (positive) eigenvalues^ Xa of A^A. The authors noted that the critical 
behaviour is expressed in terms of the (unsealed) quantities (Tk = -^ Y.b \ f°^ 
/c = 1, 3, • ■ ■. A special role is played by ai since the model has two phases 'strong' 
and 'weak' according to ai > 2 or < 2, the two phases being split by a third 
order phase transition. Recalling the asymptotic expansions (3) (4) of the DJM 
solution, precisely this order separates the 2; > and z <0 spherical free energy, and 
suggests the identification 2 ~ ai — 2, the latter having been identified in ref.[10] 
as proportional to the 'cosmological constant'. The authors recognized that the 
simplest critical point in their model is essentially the m = l point of the ordinary 



^ It seems likely that they are related in some way to the eigenvalue path integral 
representation of the DJM solutions [24] [8] . 
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unitary matrix model (6) , which is known to be true also for the m = 1 DJM 
solution [8] [25] quite apart from the Miura map. This provides some intuition for 
the coincidence. More immediately it fixes the sign of the identification between z 
and their cosmological constant to be as above. Note that this implies that z -^ +oo 
is their strong coupling regime. The multicritical points in ref. [10] are reached 
by formally setting cr^ = ■ ■ ■ = a"2m-i = 0; in addition in the continuum limit all 
(T2n+i with n > m scale to zero. The authors find that there are contributions 
to all genus in the strong phase, while in the weak phase there are contributions 
only at genus g = 1 mod m. This agrees with the DJM solution[8] on specializing 
to an m-critical point tk oc dk,m (cf.eqns.(3)(4) with F = 0). This suggests the 
identification tk ~ o'2k+i — 25fc,o up to scaling which, given the form of a2k+i 
in terms of A^, is nothing but a Kontsevich-Miwa transformation [26] [11]. This 
transformation was performed for the unitary matrix model in ref. [11] and we refer 
the reader there for details. Indeed if we set tk = 2fc+i (^2fc+i — 25fc,o) then from 
eqn.(3.7) of ref.[ll] it is clear that their Dyson-Schwinger equations are precisely 
those of the DJM solution[9]: the Virasoro constraints L^t = for n > 0. (The 
L-i constraint is tied up in a similar way to that described in the previous section, 
due this time not to open strings but to a purely non-perturbative phenomenon 
which is described further in refs.[7][22][27]. The non-perturbative parameter a 
which governs the behaviour, corresponding to an eigenvalue representation of 
the DJM solution on [a, oo), is zero in the present case (cf. footnote 6). The 
authors of ref. [11] do not in fact take the continuum limit but rather analyse 
these equations around a topological point. This is not necessary. Our equations 
of course refer to the case where the continuum limit has already been taken[8] 
but this does not alter the form of the constraints since they and the t^'s scale 
homogeneously in the continuum limit.) Finally M.Newman has compared the first 
few coefficients of the expansions for the multi-critical points of the two theories and 
found precise agreement [28]. Given all the asymptotic coincidences, the matching 
Virasoro constraints are persuasive evidence also for non-perturbative agreement, 
however to completely establish equivalence we would need to show for example 
that the solutions satisfy the same differential equation(s), together with uniqueness 
properties given boundary conditions. 
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Accepting the equivalence, Gross and Newman's model suggests some interest- 
ing interpretations for the DJM solutions. We referred above to the spherical level 
third order phase transition at z = 0. However the full non-perturbative solution 
[9] has no discontinuous transition so that it is in reality at least an infinite order 
transition, if it is correct to think of this as a phase transition at all. Similar 
comments have been made in ref. [29] . Rather it has some similarity to a 'roughening' 
transition for the strong coupling phase world-sheet expansion, a concept familiar 
in lattice strong coupling perturbation theory [30]. We may think of the point 
z = as being the point where the string fluctuations become so large that the 
concept of these strings breaks down: all world-sheet observables are singular there 
although the exact theory suffers no phase transition. One might wonder if there 
is nevertheless some sort of string description in the weak phase. Interestingly 
this question may be answered by the m = 'topological' point, a well-defined 
expansion about a free integral over U{N) identified in ref. [11]. It corresponds to 
setting all tk = except to, whence the string equation^ (1) (with F = 0) yields the 
expected trivial result u = —l/{4z^). It follows that in general the ikj^O 'topological 
expansion' about this point corresponds to the 2 ^ — oo asymptotic expansion of 
(1) (e.g. (4) withr = 0). 

Summary. 

We have shown how the string-theoretic solutions of the mKdV and KdV 
hierarchy are unified by the Miura transformation in the general picture of open- 
closed string theory, establishing connections between the non-singular solutions 
of each. It would be interesting to know whether this geometrical picture extends 
more generally to the KP and mKP hierarchy. We have also presented overwhelming 
evidence that the unitary matrix model in a tuned external field and its formal multi- 
critical generalisations matches the known non-singular closed string solutions of 
the KdV hierarchy. It seems likely that lurking behind the results of this letter is a 
more general setting in terms of quantum gravity, topological gravity and integrable 
hierarchies. 



The new topological point is simply inconsistent with 7^ = 0. 
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